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ABSTRACT 

We examine how the late divergent migration of Jupiter and Saturn may have per- 
turbed the terrestrial planets. Using a modified secular model we have identified six 
secular resonances between the frequency of Jupiter and Saturn and the four apsidal 
eigenfrequencies of the terrestrial planets (^1-4). We derive analytic upper limits on 
the eccentricity and orbital migration timescale of Jupiter and Saturn when these reso- 
nances were encountered to avoid perturbing the eccentricities of the terrestrial planets 
to values larger than the observed ones. Because of the small amplitudes of the j = 2, 3 
terrestrial eigenmodes the 52 — and — 1^5 resonances provide the strongest con- 
straints on giant planet migration. If Jupiter and Saturn migrated with eccentricities 
comparable to their present day values, smooth migration with exponential timescales 
characteristic of planetesimal-driven migration (r ~ 5-lOMyr) would have perturbed 
the eccentricities of the terrestrial planets to values greatly exceeding the observed ones. 
This excitation may be mitigated if the eccentricity of Jupiter was small during the mi- 
gration epoch, migration was very rapid (e.g., r < 0.5Myr perhaps via planet-planet 
scattering or instability-driven migration) or the observed small eccentricity amplitudes 
of the j = 2,2) terrestrial modes result from low probability cancellation of several large 
amplitude contributions. Results of orbital integrations show that very short migra- 
tion timescales (r < 0.5 Myr), characteristic of instability-driven migration, may also 
perturb the terrestrial planets' eccentricities by amounts comparable to their observed 
values. We discuss the implications of these constraints for the relative timing of ter- 
restrial planet formation, giant planet migration, and the origin of the so-called Late 
Heavy Bombardment of the Moon 3.9 it 0.1 Ga ago. We suggest that the simplest way 
to satisfy these dynamical constraints may be for the bulk of any giant planet migration 
to be complete in the first 30-lOOMyr of solar system history. 
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Introduction 



During the last two decades, observational discoveries both external to and within the solar 
system have provided strong evidence that giant planets experience large-scale orbital migration af- 
ter their formation. Outside the solar system the peculiar small orbits of the so-called 'Hot Jupiters' 
and their unlikely in situ formation suggest that th ese gas giants u nderwent large-scale, radial mi- 
gration to their observed short-period orbits (e.g., iLin et al.l 119961 ). At early tii nes l-lOMyr) 



this might be accomplished via interactions with the progenitor gas disk (e.g., iPapaloizou et al 
20071 ) . At later times migratioii may be driven by scattering and redistribut ion of a planetesimal 



disk (e.g.. [Fernandez &: Iplll984l : Murray et al.lll998l : iHahn Sz Malhotralll999l ). In addition, multi- 
planet systems may suffer orbital instabilities in which giant planets undergo mutual scattering 
events that may eject planets from the system altogether and leave the surviving ones with large 
orbital eccentricities and inclinations. Such stochastic evolution can account for the orbital charac- 



teristics of many observed extrasolar plan ets (e.g.. lRasio &: Fordlll996l : IWeidenschilling &: Marzari 



19961 : iLin Idalll997l : iLevison et al.lll998l ^ 



Closer to home, the resonant structure of the Kuiper Belt testifies to the large-scale radial 
redistribution of the solar system's giant planets more generally. The eccentric and inclined orbit 
of Pluto, and the cohort of Plutinos in Neptune's exterior 3:2 mean motion resonance, can be 



natu rally explained by the outward migration of Neptune (jMalhotral 119931 . Il995l : iLevison et al 



20081 ). As the number of solar system characteristics that can be attributed to processes associated 
with planet migration continues to increase, there is a growing consensus that large-scale planetary 
migration did occur in the solar system and that many of its dyn amical structures were determined 



durin g, and represent artifacts of, this epoch (e.g., see reviews by ILevison et al 
200?! ). 



2007 



Chiang et al 



Due to their large masses and angular momenta, the giant planets are the dominant perturbers 
of the solar system. Therefore the period of giant planet migration was one in which enormous 
orbital perturbations swept through the solar system. Clearly the timing and duration of this epoch 
informs the evolution of the rest of the system. 

Because of their hydrogen and helium rich compositions, the formation of Jupiter and Saturn 
must have been completed during the first few million years of the solar system's history whil e 
the gaseous protoplanetary disk still persisted (e.g., see the review bv lLissauer &: StevensonI 20071) . 
The accumulation of Uranus and Neptune, is considerably less constrained (iLissauer et al.lll995l) 



and remains an area of active research (e.g.. iBrvden et al.ll2000bl : IWard et al.ll2001 



Thommes et al 



20031 : IColdreich et al.l 12004 : ILevison &: Morbidellil 120071 ) . Estimates of the formation timescale of 
these planets refle ct our modest knowledge of their accuin ulation and span a wide range from 



C(106-109) years (|Goldreich et al 



2004 



Lissauer et al.lll995l ). 



Virtually all models suggested for the accretion of Uranus and Neptune are inefficient at incor- 
porating the mass of an initial planetesimal disk into the ice giants on reasonably short timescales 
(e.g., <10 Myr) and the scattering of these planetesimals may lead to radial spreading and migra- 
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tion of the growing planets ([Fernandez &: Iplll984l ). As a result, disk masses, several t imes that of 
the observed planets (i.e., the so-called minimum mass nebula: IWeidenschillingi 119771 ) and forma- 
tion locations well interior to their current orbits are often invoked to account for the accretion of 
Uranus and Neptune. 

However, local concentrations of planet building blocks at planet trapping radii (such as the ice 
line or the interface between outward and inward migration) may dec rease the formation timescale of 



cores with sufficient mass (> lOM^) to initiate rapid gas accretion (jKretke &: Linll2007l : iLvra et al 



20ld ). The formation of Jupiter leads to gap formation in the solar nebula. Near the outer edge 
of this gap, surface density and pressure maxima (at radii near or beyond Jupiter's 2:1 MMR) also 
facilitate the accumulation of grain s, protoplanetary embryos and the formation of additional giant 
planet cores ( Bryden et al.l l2000al lbl). This sequential formation process may provide a scenario 



for the formation of the gas and ice giant planets before the depletion of the gaseous nebula, with 
modest, rather than extensi ve migration, rn ostly driven by their tidal interaction with both residual 
disk gas and planetesimals (llda et al.l 120001 ). 



The epoch, mode, and characteristic timescale in which the solar system's giant planets mi- 
grated has important implications for the evolution of the planetary system. Giant planet migration 
models in particle disks can be loosely divided into two classes, based on the principal mode of 
migration. Planetesimal-driven migration results from the efficient scattering of small bodies by 
giant planets and is characterized by the smooth radial divergence of planets embedded in a plan- 
etesimal disk. Recent models of this process in the solar system suggest that the giant planets' 
orbits diverge with approximately exponential migra,tion t imescales of r ~ a/d ~ 5-20Myr (e.g., 



Hahn fc Malhotra 



1999 : 



Ida et al. 



200C 



Gomes et al.l 12004 ). 



In contrast to smooth radial planetesimal-driven migration, a global instability among the 
giant planets, involving close encounters and gravitational scattering between planets, results 
i n large stochastic iurnps in the semimajor axes, eccentricities, and inclinations of the planet s 
(|Thommes et al.lll999l . 1200.4 iTsiganis et all 1200,4 iThommes et al.lboOSbl : ISatygin fc Brown! boiol ). 
Because the planets traverse large changes in semimajor axis very quickly, the bulk of planetary 
migration may be complete in a few million years. 

A recent effort to construct a comprehensive model of planetary migration has succeeded in 
explaining several distinct dynamical characteristics of the solar system as a product of specific 
giant planet migration histories. The so-called Nice model suggests that the orbits of the giant 
plan ets evolved through phases of both planetesimal-driven migration and planet-planet scattering 
(see lLevison et al.l 120071 . for a review). In initial versions of this model, the instability of the giant 
planets is initiated by the divergent mig ration of Jupiter and Saturn across their mutual 2:1 mean 
motion resonance (jTsiganis et al.ll2005l ). Dynamics associated with the crossing of the 2:1 mean 
motion resonance (hereafter MMR) play a fundamental role in the Nice model as they are also 
responsible for facilitat ing the capture of Jupiter's Trojan asteroids with large orbital inclinations 
dMorbidelh et ahlboosl l. 
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In the inner solar system, planetary formation may have been completed prior to the time when 



tion (e.g., 


Chambers &: Wet 


leril] 


1998 




Agnor et al. 


199C 


): 


Chambers 


2001: 


Ravmond et al. 


2004: 


O'Brien et al. 2006: 


Hansen 


2009 


Ravmond et al. 


2009 




VIorishima et al. 


201C) and cosmoc 


hemi- 
n was 


cal evidence (e.g.. 


Kleine et al. 


2002; 


Yin et al. 


2002) 


suggest that terrestrial planet formatio 



completed in ~ 30-lOOMyr. 

Uncertainties in the formation timescale of Uranus and Neptune (e.g. ranging from 10^ to 10^ 
yr) combined with estimates of migration timescales (r ~ 5 — 20 Myr for planetesimal-driven 
migration or much shorter for planet-planet scattering) suggest that migration of the solar system's 
giant planets may have taken place before, during, or after terrestrial planet formation. Identifying 
when and how giant planet migration occurred remains an outstanding problem. 

Indeed, the late orbital evolution of the solar system's giant planets has even been invoked 
to account for cratering events on the Moon, an explanation that necessarily requires fully f ormed 
terrestrial planets to bear witness to the epoch of giant planet migration (e.g.. lWetherilllll975l ). The 
so-called Late Heavy Bombardment of the Moon (hereaft er LHB) refers to a perio d from ^ 3.9=fc0 .1 



Ga a go when the lunar basins with known ages formed (jTera et al.lll974l : see also iHartmann et al 



2000, for a detailed review). The formation of these basins may be attributed to a sharp increase in 
the cratering rate of the Moon several hundred million years after planet formation was complete. 
The late orbital evolution of the giant planets has been identified with this dramatic increase in the 
lunar cratering rate via a variety of dynamical mec hanisms including the sweeping of resonance s 



through the asteroid belt due to planetary migration (IStrom et al 



the late formation of Uranus and Neptune (jLevison et al.l 



2001 



2005 



Minton fc Malhotrall2009l ) 



orbital instab ilities among the giant 



plane ts resulting in orbit crossing, and planet-planet gravitational scattering (jThommes et al.lll999l . 
20021 ). or an orbital instability resulting from the divergent migr ation of Jupiter and Saturn across 
their mutual 2:1 MMR as in the Nice model (jComes et al.ll2005l ). In each of these models, objects 
on previously stable orbits are dislodged by the evolution of the giant planets and delivered to 
Earth crossing orbits. Implicit to each model is the assumption that the terrestrial planets are 
unaffected by the gross instability and/or changes in the giant planets' orbits which are responsible 
for destabilizing the LHB impacting population. 

Here we examine this assumption in detail. Specifically, we examine how the divergent migra- 
tion of Jupiter and Saturn from a more compact configuration to their current orbits may perturb 
the terrestrial planets. If the terrestrial planets witnessed an era of giant planet migration, their 
present-day orbits offer two main dynamical constraints: (1) the terrestrial planets must remain sta- 
ble during the epoch of migration and (2) the dynamical structure the terrestrial planets emerging 
from the migration epoch must be consistent with the observed terrestrial system. Since the solar 
system's giant planets migrate as a coupled system, any constraints on the migration of Jupiter and 
Saturn may also be considered as implicit constraints on the formation and migration of Uranus 
and Neptune. 
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We develop linear secular models to identify migration histories of Jupiter and Saturn that 
strongly couple them to the terrestrial planets via apsidal secular resonances. We adapt models for 
secular resonant sweeping to predict the excitation of the terrestrial planets as a function of the 
orbital state of Jupiter and Saturn at the time of resonance and their migration timescale. Finally, 
we use these models to identify migration histories that are consistent with the observed orbital 
properties of the terrestrial planets. 

This paper is organized in the following fashion. In Section [2l we describe the secular model 
and its utility for examining planetary migration. In Section [3l we discuss how passage through 
linear secular resonances excites eccentricities and how these events may be used to constrain giant 
planet migration histories. In Sections H] and [5l we compare the predictions of the linear theory 
with the results of A^-body simulations. In Section [6l we evaluate model assumptions, discuss the 
implications of our prinicipal results, and summarize our findings. 

Initial results, including the development of our secular model, the identificati on of secular 



resonances, and som e numerical simulat ions were presented in a conference talk (|Agnor &: Lin 



20071 ) . Subsequently, iBrasser et al.l (j2009l ) have also examined this problem within the context of 
the Nice model using A/^-body simulations. Our linear model predictions are in general agreement 
with, and provide a theoretical framework that accounts for, the results of their simulation results. 
We compare the two where appropriate. 



2. Secular Evolution of the Planets 
2.1. Laplace Lagrange Apsidal Precession 

The apsidal orbital precession of a system of planets can be described usi ng a set of linear first 



order differential equatioiis, who se solution is a system of eigenmodes (see, e.g.. lBrouwer &: Clemence 



19611 : iMurrav &: Dermottlll999l ). The equations of motion for the evolution of the planets' eccen- 
tricities (e) and longitude of perihelia (zu) can be constructed from the disturbing function by 
assuming that terms involving the mean longitudes are small and retaining lowest-order terms in 
the eccentricity and inclinations. For planet i, Lagrange's equations for the eccentricity vector 
{hi,ki) = ei{sin Wi, cos Wi) are 



where 



nitti ^ f mi \ a , 1 , 3GMq 



4 ^ \MQ + mJ a> c^a^ 

4 \MQ + mi) a> 3/2' ^ ; 
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and rui, Oj, and rii are a planet's mass, semimajor axis, and mean motion, respectively. The 
standard terms where a = a^/a^, a> =max(aj,a;), a< =min(aj,ai), and Laplace coefficients (6p 
have been used. Post-Newtonian corrections for relativi stic precession are includ ed in the diagonal 
matrix elements where c is the speed of light (see, e.g., I Adams &: Laughlinll2006l ). 



Solving these equations for a planetary system is an eigenvalue problem whose eigenfrequencies 
{gj) are real and represent the precession frequencies of the modes. The solutions are of the form 



hi 



^eijsm{gjt + (3j] 



(5) 



fcj = ^ eij cos{gjt + f3j 



(6) 



where the values of modal eccentricity amplitudes eij and the initial phases of the modes Pj are 
constants determined by projecting the system onto some initial state at time {t = 0). 



2.2. Eigenmode Amplitudes and the Angular Momentum Deficit 

In our use of this secular model there are two convenient and complementary formulations 
of the eigenvalue solution, each with its own advantages. The modal eccentricity amplitude (ejj) 
solution (above) clearly identifies the orbital variation of each planet i resulting from the coupling 
with the other planets via mode j. As the orbital elements of a system are its principal dynamical 
coordinates, this formulation is easily interpreted and can be compared with the results of A^- 
body simulations (e.g., using Fourier analysis to identify fundamental frequencies and amplitudes). 
However, this representation of the system's modes tends to mask the structural signature of the 
eigenvectors and the physical interpretation of a mode's amplitude. 

When considering the secular evolution of an isolated system the orbital energies and semimajor 
axes of the planets are constant and the system angular momentum is conserved. These two 
constraints may be combined to identify the so-called angular momentum deficit (AMD) as an 
important integral of a secularly evolving system. The AMD is defined as 



AMD = V miUiajil - a/ 1 - cos li), (7) 



where i runs over the planets and li is a planet's orbital inclination (see, e.g., lLaskarlll996l ). Phys- 
ically, the AMD represents the angular momentum that must be added to the system to make all 
the planets' orbits circular and coplanar. Because secular interactions do not alter the semimajor 
axes of the planets, a system's AMD quantifies the amount of angular momentum that is available 
for exchange between planets via gravitational interactions and limits the amplitude of eccentricity 
and inclination oscillations of the planets. 
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The linear secular system can also be formulated with the following transformation of variables 
and matrix elements: 



Hi = J miniafei sintUj = ^j/ij (8) 



Ki = y miTLiafei cos zui = ipiki, (9) 

and Ail {ilJi/'it'i)Aii . where the factor ij^j is the square root of a planet's circular orbital angular 
momentum (see, e.g., Brouwer Clemencel[l96ll ). 



The system's solution is again expressed as a superposition of the modes 

Hi = Y^ H,j = CjVij sm{gjt + pj) (10) 
j j 

= ^ Kij = ^ CjVij cos{gjt + /3j), (11) 
j j 

where vj = {vij} is the jth orthonormal eigenvector and Cj its modal amplitude computed from a 
set of orbital elements and masses. The system quantity | mjnja?e? = ^{Hf + Kf)"^ = ^ J2j 
corresponds to the horizontal component of the system's AMD in our secular model. The secular 
modal amplitudes {Cj} are integrals of the system and reveal the fraction of the system AMD 
partitioned into a single mode. Conversion back and forth between the two formulations is easily 
accomplished using 

e, = (12) 

As the mass distribution of a planetary system evolves (e.g., due to planetary migration, stellar 
spin down, dispersal of the gaseous nebula, etc.), the system's eigenfrequencies, eigenvectors, and 
modal amplitudes become time-dependent. However, if the dynamical tuning of the system is slow 
with respect to the eigenfrequencies {gj} and the system is not near a secular resonance (i.e., Pi — Pk 



is no t small) , the variation in mode amplitudes tends to oscillate rapidly and average out (see 



Ward 



198ll . for a complete development). For slow enough changes in the system, the modal amplitudes 
{Cj} can be considered as effectively constant. We exploit this convenient aspect of this formulation 
when exploring the consequences of secular resonant sweeping during planetary migration. 



2.3. Linear Solution of the Present-day Solar System 

We have calculated the normal mode solution for the eccentricity evolution of a six-planet 
system consisting of Jupiter and Saturn and the terrestrial planets. Uranus and Neptune have 
been omitted from this treatment for simplicity as their influence on Jupiter and Saturn is weak 
and not efficiently communicated to the inner solar system. 

The system's orthonormal eigenvectors (vj) are shown in Figure [J with their associated eigen- 
frequencies (gj)- As is the convention for the solar system, the eigenmode index (j) identifies the 
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planet that has the largest component in the eigenvector (i.e., Mercury has the largest component 
in the j = 1 mode and its evolution is strongly governed by this mode, Venus has the largest 
component in the j = 2 mode, etc.). Large components for multiple planets in a single eigenvec- 
tor indicate strong coupling between those planets. The sign of an eigenvector component (vij) 
indicates the relative apsidal orientation within the mode, with like signs indicating alignment and 
opposite signs anti- alignment. 

The eigenvectors of this six-planet system can be divided into two groups, those where the 
terrestrial planets have the largest components {j = 1-4) and those where Jupiter and Saturn have 
the dominant components (j = 5-6). This trend in the structure of the eigenvectors is indicative 
of the weak coupling between the terrestrial and giant planets in their current configuration. The 
eigenvector signatures of the terrestrial modes (j = 1-4) each have large components for at least 
two planets indicating significant coupling in their evolution. 

1.0 
0.5 
m 0.0 
I -0.5 
-1.0 

1.0 

^ 0.5 
" 0.0 
i -0.5 
-1.0 

1.0 

^ 0.5 



-1.0 

Fig. 1. — Orthonormal eigenvectors of our six-planet model including Mercury through Saturn. 
Each eigenvector is labeled with its eigenfrequency. Large components for multiple planets indicate 
strong coupling between those planets in the mode. The corrections due to nearby first-order mean 
motion resonances were included. 
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Using the VSOP82 ephemeris of iBretagnonI (jl982l ). we have computed the initial amplitudes 
and phases of this system's normal modes. The fraction of the system's AMD present in each mode 
(C|/2) is shown in Figure [2j The modal eccentricity amplitudes are listed in Table [TJ The j = 5, 6 
modes of Jupiter and Saturn contain the majority of the system's AMD0 Among the modes that 



^The neglected contributions of Uranus and Neptune (and the resulting j — 7, 8 modes) amount to an additional 
2% of the system AMD. 
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control the evolution of the terrestrial planets the j = 1 and j 



of the terrestrial subsystem's AMD (i.e., the sum S^^^C?/2). The j 



4 mode amplitudes contain ~ oo7o 
'j=i_-j — J 1 and j = 4 modes strongly 

affect the evolution of Mars and Mercury and are responsible for the large eccentricities of these 
planets. Similarly, the j = 2, 3 modes strongly affect the evolution of Venus and Earth. The 
relatively small amplitudes of these modes are manifested as smaller eccentricities of Venus and 
Earth. 

This dynamical configuration and partitioning of the AMD between the terrestrial planets is 

peculiar. The terrestrial eigenfrequencies with the largest amplitudes are both in close proximity to 

an eigenfrequency with a much smaller amplitude. For example, the j = 1 and j = 2 frequencies are 

close with gi = 5.86 yr and 52 = 7.42 yr , and the amplitude ratio is C1/C2 ~ 5. Similarly, the 

// _, 

j = 3 and j = A modes that strongly affect Earth and Mars are close in value with = 17.39 yr 
and (74 = 18.04 yr and have an amplitude ratio of C4/C3 ~ 2. 



1.0000 



0.1000 



CM 0.0100 



o 



0.0010 



0.0001 




1 2 3 4 5 6 
Mode Number] 



Fig. 2. — Fraction of the solar system's angular momentum deficit (AMD) in each eccentricity 
mode is shown. Note: AMD = Cj /2. 

The frequencies of the terrestrial modes (i.e., j = 1-4) all take values between those that 
dominate the present-day evolution of Jupiter and Saturn (i.e., 55 < (71-4 < ge). As a result the 
giant planet precession frequencies need not be grossly modified to resonate with terrestrial modes 
(e.g., (75 = gi). Meteorites and asteroids placed in or near a secular re sonance with the giant planets 
may be stron gly perturbed to eccentric orbits (IFarinella et al.lll993l ). and may even be driven into 
the Sun fe.g.. Iciadman et al.lE997l : Ichambers fc Wetherilllll998l ). 



Linear secular theory has a number of properties that make it a useful tool for exploring 
the implications of giant planet migration. As described above, the secular frequencies {gj} and 
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Mode Index j 






Planet 


i 


1 


2 


3 


4 


5 


6 


Mercury 


1 


0.1808 


0.0216 


0.0164 


0.0015 


-0.0373 


-0.0162 


Venus 


2 


0.0097 


-0.0185 


-0.0166 


-0.0017 


0.0191 


0.0094 


Earth 


3 


0.0005 


-0.0119 


0.0111 


0.0202 


-0.0017 


-0.0066 


Mars 


4 


0.0010 


-0.0248 


0.0277 


-0.0673 


0.0012 


0.0100 


Jupiter 


5 


0.0000 


0.0000 


0.0000 


0.0000 


0.0433 


0.0162 


Saturn 


6 


0.0000 


-0.0000 


0.0000 


0.0000 


-0.0357 


0.0484 



Table 1: Ecc entricity A r nplitu des (ejj) Computed in Our Six-planet Secular Model Using the 
Ephemeris of iBretagnonI (119821 ) 



eigenvectors {vj} are functions of the planet masses {rrii} and semimajor axes {aj}. This allows 
us to compute the evolution of the system's eigenfrequencies, and their sweeping through the solar 
system, for various orbital migration histories of the planets. When a system's mass distribution 
does not change, its mode amplitudes {Cj} are integrals of the system's time evolution. For the 
solar system, the planetary precession periods are ~ lO'^-lO^yr. As long as the mass distribution of 
the system evolves on longer timescales and t he system is not near resonance, the modal amplitudes 
{Cj} may be considered effectively constant (jWardlll98ll ). 



2.4. Secular Precession Near First-order Mean Motion Resonances 



In the standard Laplace-Lagrange model the direct terms of the disturbing function contain- 
ing the mean longitudes and all the indirect terms in the disturbing function may be neglected 
when averaging over long timescales. However, near commensurabilities between the mean motions 
of planets, resonant terms in the disturbing function can make significant contributions to orbit 
precession. A well-known example of this in the solar system is the so-called Great Inequality 
which results from the close proximity of Jupiter and Saturn to their mutual 5:2 MMR. Jupiter 
and Saturn's nearness to this resonance modifies the (75 and gfj frequencies by ?a -j -0.5"yT~^ and 



_j_g//yp 1 respectively (see, e.g.. iBrouwer &: Clemencdll96ll : iMurrav &: Dermottlll999l ). The lowest 



order mean motion resonances encountered during the divergent migration of Jupiter and Saturn 
are the 2:1 and 3:2. Both are first-order resonances and may strongly affect the secular evolution 
of the planets. 



Malhotra et al.l (j 19891 ) developed an approach for including the influence of first-order mean 



motion resonances in the secular model and used it to account for the evolution of the Uranian 
satellites. Their treatment of the problem includes the averaged influence of the direct terms from 
the disturbing function related to first-order (i : mean motion resonances. The 2:1 is unusual 

among the first-order resonances as it also gives rise to indirect terms in the disturbing function of 
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o f; ......................... . 

7.0 7.5 8.0 8.5 9.0 9.5 
Saturn Semimajor axis (AU) 



Fig. 3. — Comparing results of modified secular theory and A^-body simulations. Jupiter orbits 
with a semimajor axis of 5.2 AU. Saturn's semimajor axis is changed from 7 to 9.7 AU. The apsidal 
eigenfrequencies of the qq and modes were computed using standard Laplace-Lagrange theory 
(gray solid line), with the direct terms of from the disturbing function of first-order mean motion 
resonances (dash-dotted black line) and with both direct and indirect contributions from first- 
order mean motion resonances (black solid lines). Also shown are results from individual A^-body 
integrations that have been Fourier analyzed to identify the characteristic precession frequencies of 
Jupiter and Saturn (gray triangles). 

Figure [3] shows the apsidal eigenfrequencies of Jupiter and Saturn as a function of Saturn's 
semimajor axis determined using several different methods. In each case, Jupiter's semimajor axis 
was held constant at 5.2 AU, the semimajor axis of Saturn was increased from just outside the 
mutual 3:2 mean motion resonance at about 6.8 AU, across the 2:1 at 8.25 AU, to 9.7 AU. 
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To determine the eigenfrequencies of the Jupiter-Saturn system as a function of their semimajor 
axes we performed a set of individual A^-body simulations, each with a different initial semimajor 
axis for Saturn. In each simulation the planets were assigned nearly circular, nearly coplanar orbits 
(i.e., e ~ 0. 005, / = e/2), rand om phase elements and were evolved for about 11 Myr. The SyMBA 
integrator (jPuncan et al.lll998l) was used in thes e simulations and the precessional effects of general 



relativity were included (e.g.. iQuinn et al 



1991 



)• 



We determined the eccentricity amplitudes (ejj) and eigenfrequencies of the system by ana- 
lyzing the evolution of a planet's eccentricity vector (/t,:, fc,:) usin g a Frequency Modified Fourier 
Transform code (hereafter FMFT; Sidlichovskv Sz Nesvornv 1997l )FI For each simulation the two 
frequencies with the largest amplitudes in the evolution of Jupiter's eccentricity vector (/15, k^) are 
shown in Figure [3] with gray triangles. The close proximity of the measured frequencies to the 
eigenfrequencies of the Jupiter-Saturn system identifies them as the 55 and qq. Since the A^-body 
simulations include all gravitational interactions between the planets, resonant, and otherwise, 
these FMFT results serve as "measurements" of the apsidal eigenfrequencies of the actual system 
and can be used to assess the fidelity of more approximate, but analytically tractable, treatments. 

As Saturn's orbit is expanded, the values of the 55 and frequencies show two important 
characteristics. First, near both the 3:2 and 2:1 MMRs the value of the qq frequency diverges. 
The behavior of the 55 is somewhat different, showing divergence behavior only near the 2:1 mean 
motion resonance. For orbital separations outside the 2:1 resonance (i.e., the semimajor axis of 
Saturn greater than 8.25AU in this case), the two eigenfrequencies decrease as Jupiter and Saturn 
continue to diverge. 

We have used secular theory to compute the modal frequencies in three different ways. The 
standard Laplace-Lagrange treatment is shown as a gray solid line in both figures. This treatment 
contains no information regarding mean longitudes of the planets and fails to account for divergence 
in the (75^6 frequencies near the 3:2 and 2:1 mean motion resonances. 

We also computed the apsidal eigenfrequ encies of Jupiter and S aturn, correcting the Laplace- 
Lagrange treatment with the direct terms of iMalhotra et al.l (jl989l ). The value of 55 and are 
shown with a dash-dotted line. The corrections due to the direct terms account for divergence of 
the 56 frequency when approaching the 3:2 or 2:1 resonances, but fail to account for the changes in 
the (75 frequency due to the 2:1 resonance. 

Finally, we compute the apsidal eigenfrequencies using the correction of both direct and indi- 
rect terms from the first-order mean motion resonances. The contributions of the indirect terms 



associated with the 2:1 MMR may be included in the model of IMalhotra et al.l (jl989l ) by adding 
—2a6j^i to their Equation (26b) where 6ji is the kronecker delta function (R. Malhotra, private 
communication 201 1)H As the 2:1 is the only first-order resonance with indirect terms it is also 



The FMFT code used in this study was derived from that used in IChristou fc Murravl (|l999l ). 
■^A more cumbersome inclusion of the indirect terms that yields similar results may also be derived from the 
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the only resonance affected by this amended theory. 

The inclusion of the indirect terms dramatically improves the agreement between the secular 
frequencies computed with FMFT-measured values and the corrected secular theory and accounts 
for the behavior of the 55 frequency around the 2:1 MMR. This model assumes that the resonant 
argument is rapidly circulating and consequently becomes less accurate near the resonant separatrix 
where the circulation rate slows. The continuous behavior of the corrected frequency as the 2:1 
MMR is crossed, its inflection and local maxima are model artifacts of these assumptions. For 
Jupiter and Saturn more than ~ 0.1 AU outside the 2:1 MMR resonance, the results of the A^-body 
simulations and the modified secular theory are in close agreement. 

The influence of additional mean motion resonances between Jupiter and Saturn are also visible 
in FMFT measurements of 55 and g^. These appear as deviation from the smooth general trends 
indicated by secular theory. The effects of the 5:3 and 7:4 MMRs when Saturn is near 7.3 and 7.55 
AU and the 5:2 MMR as Saturn approaches its current orbit at 9.54 AU are evident. These second- 
and third-order resonances are relatively weak compared to the 3:2 and 2:1 and their influence on 
the eigenfrequencies is both smaller and localized to configurations close to exact resonance. 

These results clearly demon strate that the divergent i nigration of Jupiter and Saturn across 
their 2:1 MMR, as suggested by iTsiganis et al.l (120051 ) and I Gomes et al.l (j2005l ). results in secular 
resonant sweeping of the 55 across the values of all four of the current terrestrial eigenfrequencies. 
Thus, this migration scenario may subject the terrestrial planets to potentially destabilizing secular 
resonant perturbations. 



3. Excitation Due to Secular Resonant Crossing 



Excitation via passage through linear secula,r resonances with the giant planets has been consid - 
ered previously in several contexts (jWard et al.lll976l: iHepperiheimeiillQSOl : iNaeasawa et al.ll2005ll . 



First, we describe the test particle solutions of 



Ward et al 



(jlOTol ) and iNagasawa et al. 
adapting them to the case of resonant sweeping due to planetary migrationQ We then extend this 
model to describ e resonances with the eigenmodes of the terrestrial planets using the subsystem 
approximation of IWardI (jl98ll ). 



disturbing function for osculating elements (see, e.g., Ellis fc Murray 2000l : Murray fc Dermott 19991 ) 



^See also iGomed (1997) and the recent work of iMinton fc Malhotral (201l|) that have used similar models to 
evaluate the migration induced sweeping of the ge and sg resonances through the asteroid belt. 
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3.1. Test Particle Excitation 



Consider a test particle in a system where planetary migration causes its orbital precession 
frequency (g) to pass through a resonance with one of the system's eigenfrequencies (gj). In a 
reference frame that rotates with frequency gj, the equations of motion of /i = esin(TO — gjt — j3j) 
and k = ecos(zz7 — gjt — f3j) for the planetesimal interacting with the resonant mode are 



^ - (9 - 9j)k = fij 
dk 

+ (5 - gj)h = 0, 



where 



na \ - 



rrii 



-hi 



(13) 
(14) 

(15) 



An analytic solution for resonance passage can be obtained by introducing a new independent 
variable (p = J Agdt, where Ag = g — gj and transforming Equations 1131 and 1141 into 



dh 
d4) 



k 



Ag 



1 dk ^ 
h+ — = 0. 



The particular solution of this expression is 

h = (Ci + c) cos (j) + {Si + s) sin 
k = — (Ci + c) sin (j) + (Si + s) cos i 



with 



Ci = / ( ^ ) cos 



Si 



Ag 
Ag 



sm 



(16) 
(17) 

(18) 
(19) 

(20) 
(21) 



Defining i = as the time of exact resonance {g = gj) and setting Ag = g — gj where gj = \dgj/dt\, 
— = —^Agt"^ = —^{Ag)^/Ag, the solutions for h, k following resonance passage 



the variable 
are 



/ \ 1/2 



. 



res — A^j 



Ag 



1/2 



Ci{t') + ^) cost'2 + (^5i(t') + sint 
Slit') + cost'2 _ j^Ci(t') + ^) sint'^ 



(22) 



(23) 



where Si and Ci are the Fresnel sine and cosine integrals and we have required that h,k ^ as 
t — oo. 
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Long after resonance passage, as t — t- oo, the test particle's new free eccentricity is 



(24) 



In essence, resonance passage excites a new contribution to the free eccentricity of the body. The 
amphtude of this response is determined by the rate of passage through resonance and is a function 
of the migration rate. It is convenient to use an exponentially decaying migration rate of the form 



Aa,; 



(25) 



where Aoj = Ojj — is the distance of the planet from its final semimajor axis (oi j) and r is 
the characteristic migration timescale. We treat r as a global migration constant and use this to 
parameterize migration speed. We note that Equation (I25p can also be used to convert between 
analytic constraints, that are easily expressed as an exponential timescale (r), and migration rates 
{dai/dt) which might be more easily measured from A^-body simulations of planetary migration. 
Using the exponential migration model. 



^ \ An: An: 



(26) 



where all terms except the timescale are functions of the planetary masses and semimajor axes at 
the time of resonance. This can also be used to identify a critical migration timescale (r*) required 
to excite the eccentricity to some prescribed value eres 



1 

2^ 



/ dg dgj 



\dai 



dai 



fires \ 



(27) 



where the mass and orbital state of the perturbing planets are encapsulated in /ij (see Equation 
(jl5p ). If the test particle's free eccentricity prior to resonance passage was zero, then Equation (|24p 
is the new eccentricity after passage of the secular resonance. However, when the free eccentricity 
is not initially zero, the initial secular coordinates h^, add vectorally with those obtained from 
resonance passage /ires,^res- The t wo contributions rnay interfere constructivel y or destructively 
depending on their relative phase 7 (jWard et al.lll976l : iMinton &: Malhotrall2008l ) and the resulting 
eccentricity can be determined from geometric considerations 



6/ = (^o + e?es + 2eoeres COS 7)^/^ 



(28) 



(|Heppenheimeij 119801 ). Recognizing the phase between contributions as effectively random and 
isotropic allows Equation (|28p to be used to as a distribution function of resulting eccentricities. 
For example, given amplitudes for the initial and resonant eccentricity vectors, 7 = 0, 7r/2, and vr 
yield the maximum, median, and minimum eccentricities expected. 
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3.2. Subsystem Approximation for Terrestrial Giant Planet Interactions 

Due to their much smaller mass and modal amplitudes, the terrestrial plan ets have a ne gligible 



influence on the evolution of the giant planets. The subsystem approximation of lWardI (jl98ll ) treats 



the terrestrial planets and giant planets as separate and independently evolving secular s ystems 



and then examines their interaction. A brief summary of this model is described below (see [Ward 



198ll . for the complete development). 




In this treatment the secular evolution matrix A, whose elements are defined in Equations ([3]) 
and dH), is divided into submatricies that describe coupling within each terrestrial and giant planet 
subgroup and between the subgroups 

(29) 

Here At is the terrestrial 4x4 submatrix and Ac is the 2x2 submatrix describing Jupiter and 
Saturn. The direct influence of the giant planets on the terrestrial planets' precession frequencies 
is included through the diagonal elements of Ay (see Equation ([3|)). The coupling between the 
terrestrial and giant planet modes are included through the 4x2 and 2x4 off-diagonal submatrices 
of A and are labeled as P and P. 

The solutions for each subgroup are assumed to be independent. As before we label the solution 
to the 4x4 terrestrial subgroup (At) 

I sin(,,t + ^,) I 

with j = 1-4. Similarly, the solution for the Jupiter-Saturn subgroup (Aq) is 

cos(i/„t + ()„) J 

where and u„ are the two-planet eigenfrequencies and eigenvectors, D„ and 5„ are the amplitude 
and phase of the two-planet modes, and we retain the index labels n = 5, 6 conventionally used for 
these modes in the solar system. Because the terrestrial planets only weakly affect the precession 
of the giant planets, the eigenmodes of the two- and four-planet subgroups are nearly identical to 
those of the six-planet system (i.e., f„ = Qn, for n = 5,6) shown in Figures [T] and [21 Hereafter we 
refer to the giant planet eigenfrequencies with the symbol f„. 





3.3. The Sweeping Resonance during Migration 



Figure m shows the four eigenfrequencies of the terrestrial subgroup along with the 1^5 frequency 
as a function of Saturn's semimajor axis. As before, Jupiter is held fixed near its current orbit at 



7.0 7.5 8.0 8.5 9.0 9.5 

Saturn Semimajor Axis (AU) 



Fig. 4. — Eigenfrequencies computed using the subgroup approximation as a function of Saturn's 
semimajor axis. Each frequency is labeled with its mode index j at the right. Jupiter is located at 
5.2 AU. The j = 1-4 modes that describe the terrestrial planets are shown with gray lines. The 1^5 
frequency computed with modified secular theory secular is shown with a black solid line. Black 
diamonds show the FMFT measurements of the from a series of A^-body simulations. The 1^5 
mode, associated with Jupiter, resonates with each of the terrestrial modes during the divergent 
migration of Jupiter and Saturn. 

5.2 AU. The terrestrial frequencies are predominantly determined by their mutual interactions and 
Jupiter's proximity and are only weakly influenced by Saturn's orbital position. This is illustrated 
in Figure [Has the terrestrial frequencies (fi'1-4 shown in gray) are nearly constant as Saturn's orbital 
radius expands. 

The z^5 frequency decreases from about 12"yr~^ to about 8.5"yr~^ as Saturn's orbit expands 
from just outside the 3:2 mean motion resonance with Jupiter to about 7.8 AU. The small gaps in 
the FMFT measurement for the 1^5 frequency when Saturn is at 7.3 AU and 7.55 AU result from the 
dispersion of forcing frequencies introduced by the nearby 5:3 and 7:4 MMRs between Jupiter and 
Saturn. As Saturn's orbit expands beyond 7.8 AU toward the 2:1 with Jupiter near 8.25 AU, the 
influence of the 2:1 MMR dominates over the linear secular terms and the 2/5 frequency increases 
from 8.5"yr~^ 20" yr"^- It resonates with the (73 and 54 terrestrial frequencies when Saturn is 
near 8.17 AU and 8.20 AU, respectively. 

As Saturn's orbital radius is increased outside the 2:1, the effect of the mean motion resonance 
weakens with distance and the frequency again decreases from about 20"yr~^ to 4"yr~^ as 
Saturn's orbit expands from 8.3 AU to 9.54 AU. The resonates a second time with both the 
54 and §3 just exterior to the 2:1 MMR when Saturn is near 8.30 AU and 8.33 AU. As Saturn's 
orbit is expanded further, the 1^5 resonates with the g2 and gi when Saturn is near 8.63 AU and 
8.79 AU. In total we have identifled six individual secular resonances between the z/5 frequency and 
a terrestrial eigenmode that occurs as a result of the divergent migration, of Jupiter and Saturn, 
from just outside their mutual 3:2 mean motion resonance to their present orbits. The semimajor 
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axes of Saturn for each of these res onances is listed in Table [2l In their numerical exploration 
of this problem, iBrasser et al.l (|2009l ) parameterized the giant planet configuration at the time of 
resonance with the ratio of the orbital periods of Saturn to Jupiter [Ps/Pj). Because the value of 
the period ratio {Ps/Pj) at gj — resonances is weakly dependent on Jupiter's semimajor axis 
this ratio can be used as a convenient independent variable to approximately identify resonant 
configurations. Its value for each secular resonance is listed in Table [2j 



3.4. Secular Resonances, "Speed Traps" and Critical Migration Parameters 



Using the subgroup approximation, IWardI (jl98ll ) showed that tuning the system through a 
linear secular resonance between the terrestrial and giant planets (i.e., = gj) excites the terrestrial 
mode amplitude by an amount 



a 



j,res 



D F 



2tt 



J" I 



1/2 



where A^r^^ = gj — u^- Here Fjn is an effective width of the resonance given by 



^ in 



Vj P 



(32) 



(33) 



where vj is the transpose of Vj. Equation (j32|) is the generalization of the test particle response of 
resonance passage (Equation (p^ ) to a four-planet terrestrial subsystem. Here the mode amplitudes 
Cj,Dn take the place of modal eccentricities, and the functional dependence on the migration rate 
through Agjn is identical to the test particle case. 



The AMD exchanged with the terrestrial planets through the resonance is 



C: 



j,res 



dIt 



27rF. 



^-^ \ dai dai ' 



(34) 



The excitation of the terrestrial mode depends on three principal factors: the configuration of 
the planetary system at resonance (determined by the value of the terms in square brackets), the 
amplitude of the resonant giant planet mode (-D5), and the migration timescale (r). 



It is convenient to combine the forcing amplitude and migration timescale into a single migra- 
tion parameter 

D5 Y ( 655 \2 



5,obs 



0.0433. 



(35) 



The value of the migration parameter is proportional to the AMD exchange and can be used to 
identify the family of giant planet eccentricity states (e.g., 655) and migration timescales (r) that 
produce a certain level of excitation. In this definition, we have scaled the migration parameter so 
that f = r when Jupiter and Saturn migrate with eccentricities equal to their present-day values. 
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Using the above relations we identify the critical migration parameter for each resonance (t*^) 
capable of exciting a terrestrial mode to its observed amplitude 




In a sense, passage through the secular resonance may be used as a dynamical "speed trap" to 
constrain the eccentricity state and maximum migration timescale (or minimum migration rates) 
of Jupiter and Saturn at the time of resonance. For each secular resonance we have computed the 
giant planet orbits, resonance widths (Fjn), and critical migrations parameters (t*„) for exciting the 
resonant terrestrial mode to its observed amplitude. The characteristics describing these resonances 
are listed in Table [5J Equation (|32|) indicates that resonant excitation is dependent on the slope 
dv^/daQ. The 53 — and 54 — 1^5 secular resonances occur when Jupiter and Saturn are very near 
the 2:1 MMR and the modified secular model becomes less accurate. For these, we use values for 
the 1^5 frequency gradient interpolated from the FMFT results in Figured 

The effective width of a resonance {Fjn) is determined by the projection of the resonant 
terrestrial eigenvector (vj) onto the off-diagonal matrix P that accounts for coupling between the 
terrestrial and giant planet subgroups. As these terms scale with planet mass (see, e.g.. Equation 
dU), the effective width tends to be larger for resonances with modes where more massive planets 
have larger eigenvector components (vkj)- The resonance widths are shown in Table [2j Since 
the j = 2, 3 modes have the largest eigenvector components for Earth and Venus (see Figure [TJ 
resonances with these modes have greater effective widths than resonances with the j = 1, 4 modes. 

Also, the greater effective widths of the resonances with the j = 2, 3 modes and Equation (j32p 
indicate that for a given frequency sweeping rate at resonance {Agjn), the j = 2, 3 modes experience 
a greater excitation of their amplitudes than the j = 1, 4 modes. This susceptibility to excitation 
appears in contrast to the mode amplitude partitioning of the terrestrial planets, where C2 and C3 
are smaller than Ci and C4 (see Figure [2]). 

This general model can be used to predict the excitation of modal eccentricity amplitudes due 
to passage through a secular resonance (cjj^res) 

where obs the eccentricity amplitudes of the present-day system listed in Table [TJ Note that 
in both the test particle and subgroup approximation the resonant forcing of eccentricity scales as 
i'qj. equivalently eres oc e^^r^^'^). 

These relations and the observed state of the terrestrial planets can be used to constrain the 
giant planet migration parameter at the time of passage through a secular resonance. For example, 
if we know Jupiter's eccentricity at the time of resonance (e.g., the current value) then we can 
constrain the migration timescale of Jupiter and Saturn (and vice versa). 
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Mode 


Saturn 


Period Ratio 


Frequency 


Resonance Width 


Critical Migration Parameter 


3 n 


as(AU) 


Ps/Pj 






t;„ (Myr) 


3 5 


8.17 


1.968 


17.59 


0.040 


0.598 


4 5 


8.20 


1.979 


18.34 


-0.028 


3.609 


4 5 


8.30 


2.016 


18.30 


-0.025 


1.840 


3 5 


8.33 


2.025 


17.57 


0.036 


0.329 


2 5 


8.63 


2.138 


7.51 


-0.052 


0.050 


1 5 


8.79 


2.196 


5.87 


0.013 


0.683 



Table 2: Secular Resonances During the Divergent Migration of Jupiter and Saturn (See e.g., Fig- 
ure H]) 



Note. In this model Jupiter's semimajor axis is constant at 5.2 AU. 



The critical migration parameters i'T'jn) required to excite the terrestrial modes to their ob- 
served amplitudes are listed in Table [2j All of these migration parameters are shorter than the 
migration timescales ob served in self-consistent models of planetes imal-driven giant planet migra- 
tion (e.g., r ~5~20Myr; iHahn fc Malhotralll999l : iGomes et al.ll2004l ) by a factor ranging from a few 
to 200. These secular resonances constrain migration to occur with a smaller 655 Jovian eccentric- 
ity amplitude and/or short migration timescales to avoid driving the terrestrial eccentricities to 
large values. Further, we note that resonances between the 2^5 frequency and the §2 and (73 terres- 
trial modes, predominantly responsible for the small eccentricities of Earth and Venus, place the 
strongest constraints on the migration of Jupiter and Saturn. We compare these model predictions 
with the results of orbital integrations below. 



4. Jupiter and Saturn's Migration Near the 2:1 MMR 

As Jupiter and Saturn diverge across their mutual 2:1 MMR they perturb the terrestrial 
planets via at least five distinct dynamical events. The 53 — z/5 and the 34 — ^5 resonances are each 
encountered on both sides of the 2:1 MMR. In addition, the 2:1 MMR excites the eccentricity of 
Jupiter and Saturn on a timescale comparable to the circulation period of the resonant argument. 
Near the 2:1 MMR this is 0(10'^) yr. Because the typical secular timescales of the terrestrial planets 
are more than an order of magnitude longer (> O(IO^)), exciting the eccentricities of Jupiter and 
Saturn in this way can be thought of as rapid, non-adiabatic changes in the C5 and Cq secular 
modal amplitudes. 

To clarify the implications of individual events our analysis separates the 33 — and (74 — 1/5 
secular resonances, that occur as Jupiter and Saturn migrate near the 2:1, from the (72,1 — 1^5 
resonances, that occur for larger Jupiter-Saturn separations. We analyze the dynamics of these 
individual resonances using both A^-body simulations and secular theory and examine how they 
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act in concert. 



4.1. Initial Conditions 



To clearly illustrate the role of these events using orbital integrations it is useful to start with 
terrestrial mode amplitudes (C1-4) that are initially very small. The initial conditions for these 
simulations were prepared in the following way. Jupiter was placed on its current orbit and Saturn at 
7.2 AU. Both giant planets were given initial eccentricities of 0.01 and comparable inclinations. The 
terrestrial planets were given their current semimajor axes, small eccentricities, and inclinations 
(e = 0.01) and randomly chosen angular phases. This configuration was processed further by 
integrating it forward in time while applying an additional force to damp the eccentricities and 
inclinations of the terrestrial planets with a decay timescaie of three million years. The timescaie 
was chosen to be slow compared to the secular timescaie so that i t effectively darnped t he terrestrial 
mode amplitudes rather than the elements of a single planet ( Agnor &: Ward I2OO2I ). For these 
initial conditions the FMFT-measured eccentricity amplitudes are Cij < 0.005 for i,j = 1-4 (for 
the terrestrial planets) and 655 ~ 0.015 fo r Jupiter, about one-t hird its observed value. This value 
of 655 is comparable to that suggested in iTsiganis et al.l (120051 ) . For these initial conditions, the 
critical migration timescales needed to excite the terrestrial modes to their observed values are 
about a factor of nine longer than the migration parameter listed in Tabled] (see Equation (I35p ). 

We have performed a series of A^-body orbital integrations using this initial state. In these 
models Saturn is forced to migrate outward across the 2:1, halting its migration at 8.5AU. The 
migrat ion force applied produces exponential migration with an e-folding timescaie (r) similar to 
that of lMalhotral (jl993l ). The starting and ending points of Saturn's migration leave the 1^5 frequency 
between the 32 and frequencies and avoids depositing the system near a secular resonance (see, 
e.g.. Figured]). This allows us to isolate the effects of individual events during migration and use 
FMFT analysis of the final states to measure the ultimate eccentricity amplitudes that result. 

In this set of simulations we varied the migration timescaie from r = 10^ — > 4 x 10^ years. The 
low end of this range is near the secular timescaie and illustrates the response of the system to rapid, 
essentially non-adiabatic, evolution. The upper end encompasses migrat ion timescales suggested 
by A^-body simula tions of planetesimal-driven giant planet migration (e.g.. lHahn &: Malhotralll999l : 
Gomes et al1l2004l ). 

Figure \5\ shows an example of the coupling between the gas giant and terrestrial planets as 
Jupiter and Saturn are forced to migrate across their mutual 2:1 mean motion resonance with 
an e-folding timescaie of r = lOMyr. The top frame shows the evolution of the semimajor axes, 
pericenters, and apocenters of Jupiter and Saturn. Their divergent crossing of the 2:1 MMR occurs 
near t/r = 0.425 and the excitation of the giant planet eccentricities is readily apparent as the 
increased separation of the pericenter and apocenter curves from the semimajor axes. The lower 
panel shows these same orbital parameters for the terrestrial planets. Mars' eccentricity is excited 



-22- 



to a value similar to the observed one. As Earth and Venus have significant components in both the 
J = 3 and j = 4 modes (see, e.g., Figured]), their excitation increases eccentricities of these planets 
as well. This excitation is coincident with Jupiter and Saturn migrating across the 2:1 MMR. 
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Fig. 5. — Forced migration simulation with r = lOMyr. The semimajor axes of the Jupiter and 
Saturn and the terrestrial planets are shown with a black solid line. The pericenter and apocenter 
are shown in gray. The excitation of the terrestrial j = 3, 4 modes is evident in the increased 
eccentricity of Earth, Venus, and Mars. 



4.2. The First Crossing of the — and — Resonances 

Figure [6] shows the eccentricity excitation due to the first passage through the 53 — 1^5 and 
54 — secular resonances as a function of the migration timescale r. The theoretical predictions 
for the j = 3 and j = 4 eccentricity amplitudes of Earth and Mars computed with Equation (|37p 
are shown with dotted and dashed lines. The maximum eccentricity for planet i results when 
the contributions from these two modes are in phase, yielding a value that is simply their sum 
ei,max = + 644. This is shown with a solid line for both planets. We compare these predictions 
with the results of A^-body integrations by identifying the maximum osculating eccentricity of Mars 
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and Earth observed in each simulation prior to the time when Jupiter and Saturn crossed their 2:1 
MMR. These values are shown with filled circles. 

For short migration timescales (r < 10^ yr), the initial eccentricities (cjj ~ 0.005) dominate 
over contributions from the secular resonances. For longer migration timescales, the excitation of 
the resonances is more apparent. The A/"-body results always fall below the theoretical maximum 
and above the predicted amplitudes due to single modes. The time required to reach the maximum 
eccentricity of a two-mode system is the about half the beat period of two eigenfrequencies. For the 
j = 3 and j = ^ modes this is ~1 Myr. For r > IMyr there is sufficient time in the simulation to 
realize the maximum eccentricity prior to the Jupiter-Saturn 2:1 crossing and the values measured 
from iV-body simulations cluster tightly near the theoretical maximum. For r > 10^ yr the max- 
imum eccentricities of both planets follow the r^^^ scaling of Equation (]37p and are within 0.005 
(the initial amplitude) of the theoretical prediction. Clearly, the theory described above provides 
an accurate description of excitation during first passage through the 53 — z/5 and 34 — 1^5 secular 
resonances. 

We note that the predicted excitation of both modal eccentricity amplitudes of Mars (643 and 
644) are comparable (see Figure [6)^a)). For a given migration timescale r, the rate of frequency 
sweeping through the 53 — resonance is less (see \di>^/daQ \ in Figure H]), and the effective resonant 
width (-F35) greater, than for sweeping through the 54 — 1^5 resonance. Both factors contribute to 
the greater excitation of C3 relative to C4. Finally, the eigenvector components of Mars in the j = 3 
and J = 4 modes are comparable in magnitude (i.e., |i'4,3| ^ |^^4,4|, see Figure [I]). The net eff'ect is 
that for any migration timescale (r), the eccentricity excitation of the j = 2> component in Mars' 
eccentricity is actually greater than that of the j = 4 mode (i.e., 643 ~ l.le44, see Figure (6]^ a)). 
This appears in contrast to the observed eccentricity partitioning for Mars, where 643 = 0.020 and 
= 0.067. 



4.3. Crossing the 2:1 MMR 

The increase in the eccentricities of Jupiter and Saturn as they diverge across their mutual 
2:1 resonance might be communicated to the terrestrial planets. This change in the eccentricities 
and longitude of perihelia occurs over the circulation period of the 2:1 resonant argument (0(10'^) 
yr). As this change occurs over a timescale much shorter than the secular periods, this eccentricity 
excitation might be viewed as a sudden change in system's mode amplitudes Cj and phases 
Examination of the 2:1 MMR crossing finds that this impulsive event fairly purely excites the Cg 
mode amplitude and that t he excitation of the C5 mode, is weak (655 ~ O(10~^)) and insufficient to 



explain its observed value ( Morbidelli et al. 20091 ) Fl In general, this resetting of the giant planets' 



eccentricities should alter the entire system's mode amplitudes. However, in the observed solar 



^The conference presentation of 



Cuk 



(20071) also highlighted this aspect of the 2:1 crossing. 
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Fig. 6. — Eccentricity of Mars {i = 4) and Earth (i = 3) excited during Jupiter and Saturn's 
approach to their 2:1 mean motion resonance as a function of the migration timescale (r). The 
dotted and dashed Unes indicate the prediction of resonant excitation of the e^s and ej4 eccentricity 
amplitudes. The theoretical maximum eccentricity that might result is the sum of these two 
contributions and is shown with the solid line. The maximum osculating eccentricity observed 
in the simulations prior to 2:1 MMR crossing is shown with a filled circle. In both the linear theory 
and the A^-body simulations the forcing eccentricity of Jupiter is about one-third the observed value 
(i.e., 655 ~ 0.015). The critical migration parameters f^g and f^g suggest migration timescales of 
r = 5 and r = 30 Myr can excite the j = 3 and j = 4 modes to their observed amplitudes. 

system, the j = 6 mode weakly contributes to the eccentricity evolution of the terrestrial planets 
(see the small values of e^e in Table [T|). So, the terrestrial eccentricity amplitudes (e^j, where 

i,j = 1 4) are modified by small amounts to accommodate an increase in the j = 6 amplitude. 

On this basis, we suggest that the impulsive excitation of the Cg mode amplitude is only weakly 
communicated to the terrestrial planets and its effect on the terrestrial j = 1~4 modes is small 
relative to contributions from the gj — secular resonances. 

4.4. The Second Crossing of the (73 — and 5(4 — z/5 Resonances 

The excitation obtained from the second crossing of the 54 — and 53 — 1^5 resonances can also 
be predicted using Equation ()37p . For a given migration timescale the excitation of the j = 3, 4 
resonances encountered inside of the Jupiter-Saturn 2:1 MMR are comparable to the excitation 
due to these secular resonances outside the 2:1. 

Depending on the relative phasing between excitation from the first and second resonance 
crossings, these two contributions may interfere constructively or destructively and result in eccen- 
tricities larger or smaller than predicted for passage through a single resonance. Figure [7] shows 
examples from iV-body simulations that illustrate both constructive and destructive interference 
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due to successive passages through the 53 — and (74 — 1^5 resonances as Jupiter and Saturn cross 
their 2:1 MMR. 

In both examples the migration timescale r is sufficiently long that the first secular resonance 
crossing near t/r = 0.4 excites both j = 3 and j = 4 modes to amplitudes larger than their observed 
values. For the case shown in (a) the second passage through resonance results in a further increase 
in eccentricity amplitudes. In the second example (b), the contributions from the two resonances 
interfere destructively resulting in a final eccentricity state consistent with the observed terrestrial 
planets. This occurs despite a migration parameter greater than the critical value. 

The latter result illustrates that Jupiter and Saturn's crossing of the 2:1 MMR can temporarily 
excite the eccentricity of Mars to values of 0.2 or larger. If Mars' orbit was excited in this way it 
may have facilitated the d epletion of otherwise dyna. i nically stable and u npopulated regions exterior 
to Mars' a = 1.7-2.0 AU (|Evans TabachnikI 12002 : ISottke et allboiol ). 
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Fig. 7. — Forced migration simulations with r = 20 and 40 Myr. The evolution of the Jupiter and 
Saturn (not shown here) is similar to that in Figure [5j The black and gray lines represent the 
semimajor axes, pericenter, and apocenter distances, respectively. In (a) the contributions from 
crossing the secular resonances twice has an additive effect, resulting in a large final eccentricity 
for Mars. In (b) the excitation of Mars, Earth, and Venus due to two passages through the 53 — 1/5 
and 54 — f5 resonances partially cancel and leave Mars, Earth, and Venus in a dynamical state 
similar to the observed one. Because of the smaller Jovian eccentricity amplitude 655 ~ 0.015 the 
simulation in (a) and (b) correspond to migration parameters of 2.4 and 4.8 Myr, respectively. 

For each simulation in this set of models we conducted an additional simulation to measure the 
eccentricity amplitudes (ejj) of the terrestrial planets via Fourier analysis. The FMFT-measured 
values for Earth (633) and Mars (644) are shown in Figure [8] with filled circles as a function of the 
migration timescale r. Using linear secular theory we have computed the eccentricity amplitude 
for each planet due to secular resonance crossing both interior and exterior to the 2:1 MMR. As 
expected the net effect of the contributions from these two resonances is stochastic and depends 
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on the relative phase between them. Using Equation (j28p as a distribution function we show the 
median eccentricity expected with a black solid line. The eccentricity amplitudes of the 25th and 
75th percentiles (P25 and P75) are shown with dotted lines. Each line follows the r^/^ scaling. 
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Fig. 8. — FMFT-measured eccentricities for Earth (633) and Mars (644) are shown as a function of 
the migration timescale (r). The scatter in these measured results is due to interference between 
contributions from multiple secular resonance passages. Using secular theory and the distribution 
function of Equation 1)280 we also show the eccentricities of the 50th (solid line), 25th, and 75th 
percentiles (P25 and P75 dotted lines). For comparison the principal eccentricity amplitudes of 
Earth and Mars are 633 ~ 0.01 and 644 ~ 0.07, respectively (see Tabled]). 

Due to the stochastic influence of multiple resonance passages there is significant scatter in the 
eccentricities produced. However, about half of the FMFT-measured eccentricity amplitudes lie 
between the 25th and 75th percentile curves indicating that the results are in good accord with the 
distribution suggested by Equation (I28p and demonstrate that the linear model effectively predicts 
characteristic levels of excitation and interference. 



4.5. Interference of Multiple Contributions 

If the z^5 frequency swept through resonances with both the and 54 frequencies then the large 
amplitude of the j = 4: mode relative to the j = 3 mode requires explanation. Perhaps the relative 
amplitudes of the j = 3,4 modes resulted from the interference between multiple contributions. 
Linear secular theory and our results above suggest that the final modal eccentricity amplitudes 
eij can be calculated as the vector sum of multiple contributions (as in Equation [28]) . The — 
and 54 — resonances are close together and are each encountered twice in rapid succession as 
Jupiter and Saturn diverge across their 2:1 MMR. These resonant contributions combine with any 
primordial amplitude that resulted from the planet formation epoch. 

Studies of the gas free, late stage accretion of terrestrial planets have had difficulty simultane- 
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ously accounting for the small eccentricities of Earth and Venus and the smallish masses of Mercury 
and Mars. Further, terrestrial planet formation simulations often pro duce planets with masses 0.7- 
l.OMffl with orbital eccentricities la r ger than Earth or Venus (see, e.g.. [Chambers &: Wetherilllll998l : 



Agnor et al.l Il999l : IChamberd I2OOII : lO'Brien et al.l 120061 : [Raymond et al.l l2009l ) and suggest that a 



broad range of eccentricities might plausibly arise from accretion processes. 

Here we examine the relative probability of eccentricity amplitudes combining to yield a dynam- 
ical state comparable to the observed one for the terrestrial planets. We consider the contributions 
to each mode individually below. We compute the final value of the 633 eccentricity amplitude 
from the vector sum of an assumed initial amplitude (e.g., one resulting from formation) and two 
contributions from passage through the — 1^5 resonances on both sides of the 2:1 Jupiter-Saturn 
MMR. For a given value of the migration parameter (f), the amplitudes of resonant contributions 
are determined using Equation [37l We then choose their relative orientations at random, combine 
them and compute the resultant eccentricity vector. For each value of the migration parameter 
this procedure is repeated 10^ times and the cases where the final eccentricity amplitude is similar 
to the observed state of terrestrial planets are counted. 

When examining the results for the j = 3, 4 modes, we consider final eccentricity states with 
633 < 0.02 and 0.05 < 644 < 0.09 to be "close" to the observed dynamical state of the terrestrial 
planets. While this definition is somewhat arbitrary, the first constraint is within a factor of two of 
the small 633 amplitude, and the second is within ~ 30% of the observed 644 amplitude (see Table 
[T]for observed values). 



Constraints on the Origin of 633 

In Figure [D^a) we show the percent of cases that result in 633 < 0.02 as a function of the 
migration parameter for small (0.01), modest (0.05), and large (0.10) initial modal eccentricity 
amplitudes. Note that the critical migration parameter for the g^ — resonances interior and 
exterior to the 2:1 are f^^ ~ 0.60 and 0.33 Myr, respectively. For migration parameters larger 
than 0.7 Myr, crossing each g^ — resonance contributes a component larger than the observed 
amplitude and significant cancellation is required to achieve the observed small amplitude. For 
small initial eccentricity amplitudes (0.01), the percent of outcomes with 633 < 0.02 decreases from 
100%, to about 20% as the migration parameter increases from 0.1 to 2.0 Myr. 

When either the initial amplitude (633 = 0.05,0.10) or resonant contributions are substan- 
tially larger than the observed 633 value (e.g., for f > t^^), significant cancellation between all 
three contributions is required to produce a small final eccentricity. Figure [D^a) shows that strong 
cancellation and 633 < 0.02 can happen, but is an infrequent occurrence for large initial eccen- 
tricity amplitudes 633 = 0.05, 0.10 or migration parameters greater than 2-5 Myr. On the other 
hand, the origin of the small 633 amplitude is broadly consistent with a small primordial amplitude 
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Fig. 9. — Percent of outcomes where the 633 and 644 eccentricity amphtudes are close to those of the 
observed terrestrial planets are shown as a function of the migration parameter (Equation p5|) ). 
As Jupiter and Saturn migrate across the 2:1 MMR, the 53 — 1^5 and 54 — resonances are each 
crossed twice, once on each side of the 2:1. Contributions from these resonances are combined 
with an assumed initial amplitude and the percentages of outcomes estimated using a Monte Carlo 
technique. Curves for different assumed initial eccentricities are shown. 

(i.e., 633 ~ 0.01) and migration parameters < 1-2 Myrlfl It might also be explained by the cancella- 
tion of large resonant and initial contributions, but this becomes less likely when the contributions 
are substantially larger than the observed amplitude. 



Constraints on the Origin of 644 

As above, we compute the resultant amplitude of the j = 4 mode as being comprised of an 
assumed initial amplitude and contributions from both g/^ — resonances. The critical migration 
parameters for the resonances interior and exterior to the 2:1 are ~ 3.8 and 1.7 Myr, respectively. 
The percent of cases that produce a final eccentricity amplitude with 0.05 < 644 < 0.09 is shown in 
Figure[9l^b) as a function of the migration parameter for three different assumed initial eccentricities. 

For small initial amplitudes (644 = 0.01), this mode must be excited by the secular resonances 
to values near the observed amplitude. Resonant excitation from migration parameters < 0.2 
Myr is insufficient to perturb the j = 4 mode to the observed state (see the dashed black line in 
Figure [9)^b)). For migration parameters in the range of 0.5-1.0 Myr, about 50% of cases yield values 



®Note this upper hmit on the migration parameter f is larger than t^<^ listed in Table [5] as we are considering 
amplitudes larger than observed amplitudes of 633 as acceptably close to the observed system. 
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consistent the observed 644 amplitude. As the migration parameter increases above 1.0 Myr, the 
percent of cases consistent with the observed value of 644 decreases to about 5% at about 20 Myr. 

When the initial eccentricity amplitude is larger (e.g., 644 = 0.05), excitation from migration 
parameters f < 0.4 Myr is insufficient to alter the initial amplitude out of the acceptable range 
in more than 50% of cases (see the black dotted line in Figure [U^b)). Even with a large initial 
amplitude 644 = 0.10, cancellation between the multiple contributions produces final amplitudes 
consistent with the observed one in > 20% of cases for a large range of migration parameters 
(f ~ 0.1-2 Myr, see the gray solid line in Figure [9l^b)). 

The observed 644 amplitude is sufficiently large that it can be achieved from initial amplitudes 
in the range of 0.01-0.1 and migration parameters 0.4-2 Myr with a typical occurrence frequency 
of 20%-50% (or higher depending on specific parameter values and amplitudes). 

4.6. Implications 

The small observed 633 amplitude provides a stronger constraint on giant planet migration 
than the 644 amplitude. Specifically, the primordial, pre-migration C3 mode amplitude must be 
less than or comparable to its current value (i.e., 633 < 0.01) and the migration parameter f < 2 
Myr to account for j = 3 amplitude without requiring the low probability cancellation between 
large contributions. 

If we consider the evolution of the j = 3 and j = 4: modes as independent, we can combine 
the analysis of the individual modes to ascertain the probability of these resonance producing 
a terrestrial system with both 633 and 644 close to their observed values as the product of the 
probability of two independent events. For initial amplitudes of 633 < 0.01 and 644 < 0.05 there is 
a reasonable range of migration parameters (0.6-1 Myr) where the constraints on both 633 and 644 
can individually be satisfied in > 50% of cases. While our analysis is approximate, it suggests that 
crossing the 2:1 MMR might produce a state similar to the observed eccentricity partitioning of the 
j = 3 and j = 4 terrestrial modes in ~ 25% of cases with these initial amplitudes and migration 
parameters. 

5. The g2 — 1^5 and gi — 1/5 Secular Resonances 

For eccentricities of Jupiter and Saturn comparable to their current values as they passed 
through the g2 — 1^5 and gi — 1^5 resonances, our secular theory indicates migration parameters 
f ^0.05 Myr and 0.7 Myr will excite the terrestrial j = 2 and j = 1 modes to values larger than the 
observed ones. The order of magnitude difference in the critical parameters between these modes 
is primarily due to the larger width (-F25) of the 52 — 1^5 resonance and the smaller amplitude of 
the J = 2 mode. A migration timescale of 0.05 Myr is actually less than the period of the resonant 
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frequency (~ 0.17 Myr). For such rapid evolution it is no longer appropriate to consider the secular 
mode amplitudes {Cj} as effectively constant during migration. 

We again explore the dynamics of passage through these resonances by conducting A^-body 
simulations of forced giant planet migration with timescales that span the adiabatic secular (r ~ 1 
Myr) and non-adiabatic migration (r ~ 0.01 Myr) regimes. Migration timescales as short as 
r ~ 0.01 Myr are not likely to result from smooth migration driven by the scattering of small bodies 
by Jupiter and Saturn. Such rapid migration might result from scattering with other planets or ice 
giants. In this sense our use of smooth migration with very short migration timescales may serve 
as a simple proxy for the rapid non-adiabatic evolution driven by giant planet scattering. 

5.1. Initial Conditions 

We again isolate the influence of secular resonances by constructing initial conditions where 
the terrestrial mode amplitudes are very small, the giant planets' modal amplitudes are near their 
present values, Jupiter's semimajor axis is fixed at 5.2 AU, and Saturn's initial semimajor axis of 8.5 
AU places the z/5 frequency between the 53 and g2 frequencies. We prepared this initial condition 
using a series of orbital integrations. First, we performed an integration starting with only Jupiter 
and Saturn with their observed orbits and slowly forced Saturn to migrate inward to a semimajor 
axis of 8.5AU. Because the forced migration is slow relative to the secular timescale, the D^ q modal 
amplitudes are effectively constant during migration to the new orbits. The terrestrial planets 
were then added to the system and we conducted a second integration applying slow eccentricity 
damping to decrease the terrestrial mode amplitudes. The resulting state of the Jupiter-Saturn- 
terrestrial planet system was then used as the initial condition in the numerical experiments that 
follow. FMFT analysis of this initial condition finds the terrestrial eccentricity amplitudes with 
Cij < 0.005 and the 655 amplitude of Jupiter very near its current value (i.e., 655 ~ 0.0435). Thus, 
for this set of simulations f = r. 

5.2. Results 

Figure [10] shows the eccentricity evolution of the terrestrial planets as Saturn migrates from 
8.5AU — ?• 9.5 AU with an exponential timescale of r = 2 Myr. The 52 — and gi — 1^5 resonances 
are crossed at about t/r =0.3 and 0.6 when Saturn is near 8.7 and 8.9 AU, respectively. Excitation 
of both terrestrial eigenmodes is evident in the relative timing of the maximum and minimum 
eccentricity observed in the orbits of Mercury, Venus, and Earth. Mercury's osculating eccentricity 
reaches a maximum when contributions from both the j = 1 and j = 2 modes are in phase. For this 
same orientation between modes, the eigenvector signatures for Earth and Venus are anti- aligned 
(see Figure[T]) resulting in cancellation between modal contributions and the osculating eccentricities 
of Earth and Venus simultaneously reach a minima. For this cancellation to yield eccentricities near 
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zero, the eccentricity amplitudes from each mode must be comparable. In agreement with secular 
theory, the eccentricity amplitudes of Mercury, Venus, and Earth are each driven to values in excess 
of the present-day values. Mars is only weakly coupled to the other terrestrial planets through the 
j = 1,2 modes (see Figured]) and is minimally excited via these resonances. 

Following each migration simulation we conducted a separate orbital integration and performed 
FMFT analysis of the resulting dynamical state. The modal eccentricity amplitudes (ejj) for 
Mercury, Venus, and Earth are shown in Figure [TT] as a function of the migration timescale (r). 
The predicted excitation of the j = 1, 2 modes are shown as solid and dashed lines, respectively. 
Also the FMFT-measured eccentricity amplitudes for the j = 1, 2 modes are shown with diamonds 
and filled circles, respectively. 

For migration timescales t = 0.1 0.4 Myr, the eccentricity amplitudes closely match the 
predictions of secular theory and follow the r^/^ scaling. For Mercury the excitation of en ~ 0.20 
requires r = 0.6 Myr. Note the j = 2 component of Mercury's eccentricity is excited to half that 
of the j = 1 mode (ei2 ~ O.Sen) for any timescale r. For slower migration. Mercury may be 
excited to very large eccentricity due to the combined effect of both resonances. In these cases 
general relativity and higher or der effects co nspire to increase Mercury's precession rate and the 



value of the gi eigenfrequency (iLaskarl 120081 ). For r = 0.3 — > 1.0 Myr, the FMFT measurement 
of the gi eigenfrequency increases from 5.84"yr~^ — > 6.3"yr~^- This tuning of the gi toward the 
(72 = 7.43"/yr alters the structure of both eigenvectors. As two eigenfrequencies converge, their 



eigen vector signatures become increasingly similar (see detailed discussions of this issue in IWard 



1981 



). As the value of gi increases toward the g2, the magnitude of the components of Earth 
and Venus in the j = 1 eigenvector increase and those of Mercury slightly decrease. The j = 2 
eigenvector evolves in a complementary fashion. The component of Mercury in the j = 2 mode 
(f 12) increases while the components of Venus and Earth (^22 and U32) decrease. In this case, the 
modest discrepancies in FMFT-measured eccentricity amplitude and secular theory in Figure [11] 
result from changes in the eigenfrequencies and eigenvector signatures. We note that these effects 
become important at eccentricities well above the observed values in the present solar system. 

For migration timescales r < 0.1 Myr, the entire epoch of migration may occur in less than a 
precession cycle. Such migration is not slow with respect to the z^s frequency. Since the eccentricity 
amplitudes of the terrestrial planets were all < 0.005 prior to migration. Figure [TT] shows that 
rapid non-adiabatic giant planet migration can also contribute to the eccentricity amplitudes and 
AMD in the terrestrial modes. For the j = 2 mode, the eccentricity amplitudes of ei2 = 0.052, 
622 = 0.033, and 632 = 0.025 for Mercury, Venus, and Earth are systematically achieved and appear 
independent of the migration timescale. These values correlate with the planet's distance from the 
f5 resonance in the present-day solar system and are likely related to the forced eccentricity in the 
region. These values are also greater than the observed values for the j = 2 eccentricity amplitudes 
of the terrestrial planets. 



For the fastest migration timescale explored above, this change in eccentricity amplitudes re- 
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Fig. 10. — Forced migration simulation with r = 2.0 Myr. The semimajor axes of the Jupiter and 
Saturn and the terrestrial planets are shown with a black solid line. The pericenter and apocenter 
are shown in gray (colored in online version). 
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Fig. 11. — Final modal eccentricity amplitudes (ejj) of the terrestrial planets as a function of the 
giant planet migration timescale (r). The symbols indicate FMFT measurements of the amplitude 
from A'^-body simulations and the lines indicate the prediction of secular theory for the excitation 
of the j = 1, 2 modes in each planet (Equation (|37p ). Results for the j = 1 mode are shown with 
solid lines and open diamonds and results for j = 2 mode are shown with dashed lines and filled 
circles. 



suits from a single rapid jump from one orbital state of Jupiter and Saturn to another. However, 
giant planet migration due to scattering between the planets would involve se veral to many stochas- 
tic jumps in semimajor axes and modal amplitudes (jMorbidelli et al.l 12003 ). Repeated scattering 
events among the giant planets may induce stochastic and diffusive exchange of AMD between the 
terrestrial and giant planet mode amplitudes. If the eccentricities shown in Figure [TT] are indicative 
of the size of the steps, then this type of evolution may also perturb the terrestrial system to ec- 
centricities larger than the observed values. As in the case of slower planetesimal-driven migration, 
selective cancellation between multiple contributions may be required for the terrestrial planets to 
emerge from this type of giant planet evolution with an ecc entricity state resembling the observed 
one, even if the initial j = 2, 3 mode amplitudes are small. iBrasser et al.l (|2009l ) show an iV-body 
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simulation result that illustrates this type of evolution. A systematic exploration of the diffusion of 
AMD between the giant and terrestrial planets during an epoch of giant planet scattering is needed 
to more generally evaluate the influence this style of giant planet migration has on the terrestrial 
planets. 



5.3. Interference of Resonant and Primordial Amplitudes 

There is about an order of magnitude difference between the critical migration parameters 
for the j = 2 and 1 modes (f^g = 0.05 Myr and f^g = 0.7 Myr, see Table [5]) with both values 
shorter than migration timescales suggested for planetesimal-driven migration by at least an order 
of magnitude. Could cancellation between the resonant contributions yield a dynamical state 
consistent with the terrestrial planets for migration timescales of 2-10 Myr? To evaluate this 
hypothesis we model the addition of multiple contributions to the j = 2, 1 modes using the Monte 
Carlo technique described in Section 14.51 Because the §2 — 1^5 and 171 — z/5 resonances are crossed 
once during the last ~ 1 AU of divergence between Jupiter and Saturn, the final mode amplitude 
results from just two contributions, the secular resonant one and the primordial amplitude. When 
comparing results with the terrestrial planets, we consider "close" to the observed dynamical state 
to mean 622 < 0.03 and 0.12 < en < 0.23 for the j = 2 and j = 1 modes, respectively. As above, 
these ranges are somewhat arbitrary, but are comparable to the observed 622 = 0.02 amplitude and 
within 30% of the observed en = 0.18 amplitude. 
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Fig. 12. — Percent of outcomes where the e22 and en eccentricity amplitudes are close to the 
observed dynamical state of the terrestrial planets are shown as a function of the migration param- 
eter (Equation (j35|) ). Contributions from the g2 — and gi — 1^5 resonances are combined with 
an assumed initial amplitude as in Figure O Curves for different assumed initial eccentricities are 
shown. 



- 35 - 



Figure [T2] shows the percent of outcomes consistent with the observed 622 and en as a function 
of the migration parameter for several different assumed primordial eccentricity amplitudes. Can- 
cellation can produce systems with 622 < 0.03 for migration parameters larger than the critical one 
(see, e.g., Figure [T2](a)), but the requisite destructive interference requires primordial and resonant 
contributions of comparable size. Also, the probability of the required cancellation decreases as 
the size of contributions (and migration parameter) increases. For example, migration parameters 
of 1 Myr contribute a vector of about 0.10 in magnitude to 622 (see Equation (1371 and Table [2]). 
Combining this with an initial eccentricity amplitude of 0.10 will produce a sufficiently small final 
value of 622 in ^ 10% of cases. 

Figure [T2rb) shows results for the en amplitude. Cancellation with an initial en amplitude 
generally broadens the range of migration parameters that are able to produce systems close to 
the observed state of the terrestrial planets. However, even large initial eccentricities of 0.10 do 
not permit migration parameters much larger than 1 Myr, and are not consistent with migration 
timescales suggested for planetesimal migration, unless the amplitude of the ess amplitude is lower 
than its observed value as these resonances are encountered. 

Recall that the iV-body simulations with migration timescales r = 0.01 Myr are fast with 
respect to the secular timescale. Also for very fast migration the resulting en and e22 amplitudes are 
comparable to those obtained by migration parameters f ~ 0.1 Myr. We use this level of excitation 
to estimate the minimum forcing expected during instability-driven migration (i.e., planet-planet 
scattering). For example. Figure [T2l shows that for modest initial eccentricity amplitudes (i.e., en ~ 
e22 — 0.05) and fortuitous interference between primordial and resonant contributions, both the 
large observed value of en and the small observed value of e22 might be produced in ~ 10% of cases 
for migration parameters of 0.10-0.4 Myr. This fast exponential migration modeled here might be 
thought of as the rapid divergence of Jupiter and Saturn by a single strong scattering event, or 
one large jump in semimajor axis to the observed orbit. When Jupiter and Saturn experience 
multiple scattering events the forcing of the terrestrial planets may be greater and the success rate 
in producing the observed state of the inner planets lower. 

Because the en modal amplitude is about an order of magnitude larger than e22 and both 
are excited in comparable ways by the sweeping resonance (see Equation (f37|) , Table [21 and 
Figure [TT]) , it is difficult to explain their very different amplitudes as a natural outcome of plane- 
tary migration, without invoking strong interference between multiple contributions or additional 
processes. 



6. Discussion 



Using a modified secular theory we have calculated how the i^s and eigenfrequencies of the 
Jupiter-Saturn system evolve as these planets diverge from just outside their mutua l 3:2 m ean mo- 
tion resonance to their observed orbits. Using the subgroup approximation of I War dl (Il98ll ) we have 
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identified six secular resonances between the frequency and one of the terrestrial eigenfrequencies 
during this divergent migration. 

Based on analytic arguments we combined the Jovian eccentricity amplitude (655) and orbital 
migration timescale (r) into a single migration parameter f = (e55/0.0433)^r that encapsulates the 
forcing strength experienced by the terrestrial planets via secular resonances with the fs frequency. 
We used this model to identify the critical values of the migration parameter required to excite 
the terrestrial eccentricity eigenmodes to their observed amplitudes. We validated the analytic 
secular model with A^-body simulations and showed that linear secular models can account for 
the fundamental character and scaling of the eccentricity excitation experienced by the terrestrial 
planets via these secular resonances. We then used this model and the observed modal amplitudes 
of the terrestrial planets to constrain the migration of Jupiter and Saturn when these secular 
resonances were encountered. 



6.1. Evaluating Model Assumptions 

To clarify the dominant processes at work, we have made a number of simplifying assumptions. 
First, we have neglected the influence of Uranus and Neptune. The distant locations and relatively 
small mass of the ice giants render them weak perturbers of Jupiter and Saturn. Their omission 
from the secular model does not significantly modify the 1/5 and vq frequencies of Jupiter and Saturn 
or those of the terrestrial planet subsystem (e.g., Agi/gi ^ 0(10""^)). 

Similarly, we have neglected any modification of the system's eigenfrequencies that results 
from interaction with the planetesimal disk fueling migration. The modification of the z^s and uq 
frequencies by the planetesimal disk is a function of the disk mass and location. When considering 
the late migration of Jupiter and Saturn, the planetesimal disk must have survived for at least the 
30-lOOMyr required to form the terrestrial planets . Bodies with orbits near Jup i ter and Saturn have 



a me dian dynamical lifetime of about 6-9 Myr (jTiscareno &: Malhotral l2003l : iBailev &: Malhotra 



2OO9I ) . This suggests that the Jupiter-Saturn region was la rgely devoid of plan etesimals by the time 



the terrestrial planets were completely assembled (see, e.g.. lGomes et al.ll2005l ). A planetesimal disk 
exterior to Jupiter and Saturn with mass comparable to that of Uranus and Neptune (i.e., ~ 30M^) 
affects the 1^5 and frequencies by an amount comparable to that of the ice giant planets and can 
be similarly neglected. 

For rapid migration (r < 1-2 Myr), a large fraction of the planetesimal disk must be on 
crossing orbits with Jupiter and Saturn. To roughly estimate the effect of this disk on the giant 
planet eigenfrequencies we treat it as a single lOM^ planet between the gas giants. This increases 
both the z^s and uq frequencies by ~ 20%. However, this contribution decays as migration proceeds 
and the planetesimal disk is depleted. We suggest that these deviations in the 1/5 frequency are not 
likely to prevent resonances with the terrestrial modes or grossly modify the factors that determine 
the excitation experienced by the inner planets. 
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In exploring the problem we held Jupiter's semimajor axis constant to reduce the number 
of model parameters. More generally, the configuration of the planetary system when secular 
resonances occur depends on the semimajor axes of all the planets. So, the examples shown and 
values listed in Table [2] are not unique. Self-consistent A'^-body simulations of planetesimal-driven 
migration of the solar system typically show Jupiter to migrate inward a few tenths of an AU 



and Saturn migra ting outward a few AU on the same timescale (e.g., iHahn &: Malhotra 



199S 



Gomes et al.l 120041 ). 



We have used the secular analytic model to examine how Jupiter's migration affects the plane- 
tary configuration where resonances occur and the critical migration parameter for each resonance. 
Considering Jupiter with an initial semimajor axis of 5.4 AU, Saturn starting just outside the 3:2 
MMR with Jupiter, and both planets migrating to their observed orbits, all six of the secular res- 
onances reported in Table [2] are encountered, albeit with slightly different values of the planetary 
semimajor axes, resonant frequencies, and critical migration timescales. While the 1^5 g frequencies 
are strongly affected by the orbital divergence of Jupiter and Saturn, the effect of Jupiter's inward 
migration of ~0.2 AU on the terrestrial eigenfrequencies ((71-4) is small and only slightly larger 
than the effect of Saturn's outward migration (see, e.g.. Figured]). Consequently, the orbital pe- 
riod ratio between Saturn and Jupiter (Ps/Pj) at a particular gj — secular resonances is nearly 
independent of Jupiter's semimajor axis, or whether it is migrating. This allows th e period ratio 



to be used as a single independent variable to identify secular resonances (see, e.g., iBrasser et al 



2009 ) and for comparison with scenarios where both gas giants are migrating. 



Further, the secular model shows that the frequency gradient at resonance due to Jupiter 
{\d{gj — i'^)/da^\) is comparable to that of Saturn {\d{gj — v^)/da^\). In this situation, the planet 
with the greater migration rate (i.e., da/dt) most strongly determines the level of terrestrial excita- 
tion and the critical migration parameter (see Equation ([36]) ). Due to its larger mass Jupiter tends 
to be less mobile and its planetesimal-driven migration rate generically slower than Saturn's. Con- 
sequently, Saturn's migration rate most strongly determines excitation of the inner solar system. 
Finally, for the cases where Jupiter migrates inward a few tenths of an AU, the critical migration 
parameters vary from those listed in Table [2] by 20%-30%. 

In total, the various assumptions we have made to simplify the model and its presentation 
modestly affect the secular system. We consider the critical migration parameters listed in Table [2] 
as broadly representative of many plausible migration scenarios suggested for the solar system. 



6.2. Constraints on the Migration of Jupiter and Saturn 

As Jupiter and Saturn diverge across their mutual 2:1 MMR they encounter the (73 — 1^5 and 
gA — 1^5 secular resonances on both sides of the resonance (i.e., for Ps/Pj — 1.97-2.03). Due to 
several factors, including the small amplitude of the j = 3 terrestrial mode, the g^ — resonances 
near the 2:1 MMR of Jupiter and Saturn constrain the giant planet migration parameter to satisfy 
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f < = 0.60 Myr. Following the 2:1 crossing, the g2 — 1^5 and gi — resonances are encountered 
as Jupiter and Saturn diverge outside their 2:1 MMR (when Pg/Pj ~2. 1-2.2, see Table [2|). In a 
similar manner the migration history considered must satisfy f < = 0.05 Myr when Ps/Pj — 2.1 
and f < fj*5 = 0.68 Myr when Ps/Pj — 2.2. We note that greater values of the migration 
parameters may also produce acceptable states for secular resonances with the terrestrial modes 
(j = 1-4). However, achieving agreement with the observed terrestrial system relies on cancellation 
of multiple contributions (see, e.g., Figures [3 and [8]) . When the magnitudes of all contributions 
are known, we have shown that the probability of producing a particular eccentricity state can be 
estimated using analytical arguments (see Equation (|28|) , Figures and I12|) . We have also shown 
that the probability of yielding the required cancellation between multiple contributions decreases 
significantly when the critical migration parameter t*^ is exceeded by a factor of two to three. With 
these caveats we suggest that these values represent characteristic constraints on the migration of 
Jupiter and Saturn due to passage through these secular resonances. Further, because the solar 
system's giant planets migrate as a coupled system, constraints on the migration of Jupiter and 
Saturn should be considered as indirect constraints on the formation and migration of Uranus and 
Neptune and the evolution of the giant planet system as a whole. 

For each gj — secular resonance identified, planetesimal-driven migration of the giant planets 
(i.e., r ~ 5-20 Myr) with the Jovian eccentricity comparable to the present value would excite the 
terrestrial system to eccentricities and modal amplitudes well beyond their observed values. The 
constraints on giant planet migration imposed by these secular resonances and the dynamical 
structure of the terrestrial planets may be satisfied in several ways. 

First, if the giant planets migrated with eccentricity amplitudes comparable to their present day 
values the strong forcing of the terrestrial planets via the secular resonances may have interfered 
destructively with an initial amplitude of the j = 1-4 modes (e.g., from the accretion epoch). 
This could deposit the terrestrial planets in a dynamical state similar to the observed one. The 
general scenario requires the cancellati on of two large coni ponents and we have shown that this is 



accordingly of low probability (see also iBrasser et al.l l2009l ) . 



Using the resonant excitation model and geometric arguments we can identify additional con- 
straints required on this type of evolution. For cancellation to be viable, the initial and resonant 
contributions must be of a similar size and the excitation from resonance crossing necessarily less 
than an amount that would render the terrestrial subsystem unstable. For the g^ — and 54 — 
secular resonances encountered near the 2:1 MMR, the terrestrial planets are driven to crossing 
orbits if the migration timescale of Jupiter and Saturn is r > lOMyr and the 655 has the present-day 
value. For the 52 — 1^5 and gi — v^ resonances, f > 2Myr leads to crossing orbits and instabilities (see, 
e.g.. Figure flOl) . This latter constraint indicates that planetesimal-driven migration through the 
92,1 — 1^5 secular resonances is not viable unless the Jovian eccentricity amplitude 655 is significantly 
smaller than the observed value. 



Second, if the Jovian eccentricity (655) was smaller during migration, then longer timescales 
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might be permissible. For r ~ 10 Myr, the 53 — 1^5 secular resonance constrains the Jovian ec- 
centricity amplitude 655 < 0.011, or about a quarter of the current value, when Jupiter and Sat- 
urn migrated across the 2:1 MMR. This eccentricity is comparable to those typical of models of 
planetesimal-driven migration. However, since it is less than the observed 655 = 0.0433 value it 



also re quires that the -D5 mode amplitude be excited later via some other process. iMorbidelli et al 



(j2009l ) examined the origin of the mode amplitude in detail and found that repeated planet- 
planet scattering between Jupiter or Saturn and an ice giant planet may explain the observed 
amplitude and account for its excitation via an instability in the giant planets' orbits after the 
2:1 MMR of Jupiter and Saturn i s crossed. This analytic result is consistent with the numerical 
simulations of lBrasser et al.l (j2009l ) and the conclusions drawn from them. Using the secular model 
we have also shown for initial eccentricity amplitudes 633 < 0.01 and 644 < 0.05 the eccentricity 
partitioning of the j = 3 and j = 4 terrestrial modes can be expected with a roughly 25% prob- 
ability for migration parameters in the range f =0.6-1.0 Myr. Again, these values are broadly 
consistent with those of planetesimal-driven migration. 

When examining the (72 — 1^5 and gi — resonances, a Jovian modal eccentricity amplitude of 
655 < 0.003 is required to allow a planetesimal-driven migration timescale of r ~ lOMyr through 
these resonances without exciting the j = 2 terrestrial mode amplitude beyond its observed value. 
This value of 655 is quite small, and again requires the late excitation of the Jovian eccentricity 
(e.g., via scattering with a > lOM^ body to excite the -D5 amplitude dur ing the last 0.5-0.8 AU of 
divergent migration between Jupiter and Saturn; IMorbidelli et al.l 12003 ). 



Third, rapid, stochastic migration in large steps in semimajor axis is also possible when the 
system of planets is unstable and planet-planet scattering ensues. Such global instabilities may 
allow secular resonances to be jumped over and the whole migration epoch to be much shorter 



than suggested by planetesimal-driven migrat ion. (jThommes et al. 



Morbidelli et al. 



2007 
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Tsiganis et al. 



2005 



Batvgin fc Brownl l20ld ) . If such rapid migration was driven by scatter- 



ing between planets, this process must cause the orbits of Jupiter and Saturn to diverge. For 
Saturn to move outward it must scatter an ice giant inward. Similarly for Jupiter to move inward, 
it must scatter an ice giant outward. This style of instability-driven migration might allow the 1^5 
and secular r esonances to move q uickly through the inner solar system. A^-body simulations of 
this process by (jBrasser et al.ll2009l ) suggest that it is possible that the terrestrial planets might 
emerge from this style of planet migration with orbits comparable to their observed ones. 

On the other hand, our simulations with very fast migration (r < 0.2 Myr) indicate that the 
j = 2 terrestrial mode can be excited to amplitudes in excess o f the observed values. If the D5 
mode amplitude originated via repeated planet-planet scattering (jMorbidelli et al.ll2009l ) it appears 
plausible that instability-driven giant planet migration may perturb the terrestrial planets to an 
overly excited state. In this regard, additional work assessing the relative frequency of characteristic 
outcomes (e.g., those similar to the observed state) would be welcome. 



Finally, if planetary migration was rapid as suggested by the constraints of these secular reso- 
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nances (r <5-10 Myr), then entire epoch of giant planet formation and migration, could have been 
largely complete before terrestrial planet accretion was finished. In this case the dynamical structure 
of the inner solar system was determined later (e.g., by processes related to accretion and evolution) 
and the formation timescale of the terrestrial planets (i.e., 30-100 Myr) cons trains the timing by 
which giant planet formation and migration was largely complete. Recently, IWalsh Sz Morbidelli 
(|201ll ) have investigated the effect of the early planetesimal-driven migration of Jupiter and Saturn 
on terrestrial planet formation. Their A^-body accretion simulations with a migration timescale of 
r = 5 Myr for Jupiter and Saturn produce terrestrial planet systems with AMDs comparable to the 
observed terrestrial planets. While several open problems of terrestrial planet formation persist in 
these simulations (e.g., explaining the small mass of Mars and the clearing and dynamical structure 
of the asteroid belt), these results suggest that effects of the sweeping resonance through the 
terrestrial region may be tempered by terrestrial accretion dynamics. 



6.3. Implications for the LHB 



In the Nice model explanation of the LHB (jGomes et al.ll2005l ). an instability among the giant 
planets is initiated by Jupiter and Saturn crossing their mutual 2:1 MMR. This instability may be 
delayed by the very slow divergent migration of Jupiter and Saturn of a few tenths of an AU toward 
the 2:1 MMR over 600-700 Myr. This slow migration roughly corresponds to exponential migration 
timescales of r > 3 x lO^yr. The secular resonant model presented here predicts that crossing the 
— 1/5 and 54 — 1^5 resonances this slowly requires the Jovian eccentricity amplitude 655 < O(10~^) 
to leave the terrestrial planets in a state comparable to the observed one. A modestly larger Jovian 
eccentricity (655 > 0.001) may perturb the terrestrial planets to crossing orbits and drive a global 
instability of the terrestrial system. 

Requiring such small values of the Jovian eccentricity is a very strong constraint on the evo- 
lution of the gas giants and we interpret this as an indication that the slow divergent migration 
of Jupiter and Saturn toward the 2:1 MMR is not generally consistent with the coexistence and 
dynamical s tructure of the terr estrial planets. As a result, th e spec ific LHB initiating scenario 
described in iGomes et al.l (|2005l ) appears unlikely. iBrasser et al.l (|2009l ) have drawn similar conclu- 
sions based on the results of A^-body simulations. 

If the LHB resulted from the large scale migration of the giant planets, then an instability and 
giant planet scat tering may be requ ired to avoid strongly perturbing the terrestrial planets via sec- 
ular resonances (jBrasser et al.ll2009l ). Alternate modes of triggering a late instability (e.g., Uranus 
and Neptune's divergent migration across a m ean motion resonance) are possible and may offer 
viable alternative modes of initiating the LHB ([Morbidelli et al.l 120071 ). 



However, the j = 5 mode of Jupiter and Saturn has components in terrestrial planets com- 
parable to those of j = 2, 3 terrestrial modes (see, e.g.. Table [I]). If giant planet scattering is 
responsible for exciting the j = 5 mode amplitude to the observed value, then this excitation may 
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also be communicated to the terrestrial planets (e.g., stochastic diffusion of AMD between j = 5 
mode and the terrestrial j = 1 — 4 modes). Our simulations of very fast migration (i.e., r < 0.10 
Myr) suggest that excitation of the terrestrial planets via this process is likely, but more work 
to examine the net exchange of AMD between the terrestrial and giant planet systems is needed 
to assess the dynamical implications of instability-driven giant planet migration for the terrestrial 
planets. 



6.4. Summary 

We have shown that if the terrestrial planets witnessed an epoch when the gross orbital struc- 
ture of the giant planets changed, then the dynamical structure of the inner solar system was likely 
altered. In this scenario, the ultimate dynamical state of the terrestrial planets is a product of 
several processes (e.g., accretion dynamics, passage through secular resonances, AMD diffusion via 
giant-planet scattering, and chaotic diffusion of the system's AMD over long timescales) that act 
in concert and whose individual influences are challenging to isolate and disentangle. 

The small amplitudes of the j = 2, 3 terrestrial eigenmodes, that are closely associated with 
the small eccentricities of Earth and Venus, provide strong constraints on the migration of Jupiter 
and Saturn. While more work examining the characteristic forcing of the terrestrial planets by 
instability-driven migration is needed, explaining the observed dynamical state of the terrestrial 
planets as a natural and frequent outcome of the late migration of the giant planets remains a 
challenge. 

An alternate and perhaps simpler possibility is that giant planet formation and migration 
was largely complete before the dynamical structure of the inner solar system was determined. 
Formation models and isotopic evidence suggest that terrestrial planet accretion was complete 
in 30-100 Myr. The bulk of giant planet migration could easily be accommodated in shorter 
time intervals via some combination of tidal interaction with a gas disk, planet-planet scattering, 
and planetesimal-driven migration. Any strong perturbation of the inner solar system resulting 
from giant planet migration might then be tempered through a variety of processes that depend 
on the state of the forming planets and the protoplanetary disk at the time of the disturbance. 
Additional mechanisms that may act at this earlier time include dynamical friction with a disk of 
planetesimals, orbita l damping due to interactions with a remnari t gas disk, or a combination of 



dynamical processes (jNagasawa et al.ll2005l : iThommes et al.ll2008al ). 



In the sequential accretion scenario, the formation timescale for giant planet cores (~ lOM^) 
in the outer regions of a minimum mass nebula appear to be considerably longer than the obser- 
vationally inferred depletion timescale of protostellar disks. Ice giant formation at smaller orbital 
radii (and shorter orbital periods) is often invoked to overcome this problem, but requires the 
subsequent outward migration of the ice giant planets. However, it is also possible that disk gaps 
near Jupiter and Saturn's orbits may have provided migration barriers which eventually led to the 
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accumulation of planet-building materials and the emergence of Uranus and Neptune within a few 
million years (see Section [T|) . 

Th e resonant struc t ure of the Kuiper Belt is often interpret ed as a product of planetesimal- 



(Levison et al 



driven (|Malhqtral 1 19951 : llda et al.l l200Cl : iHahn &: Malhotral l2005l ) or instability-driven migration 



20081 ) of the giant planets. However, the giant planets' outward r nigration may also 



be in duced by their tidal interaction with a viscously expanding gaseous nebula (jLin &: Papaloizou 



19861 ) or by a disk which is undergoing photoevaporation. In these cases, giant planet migration 
may be mostly complete in '--'10 Myr and could predate the final assembly of the terrestrial planets. 
Such an early migration scenario would avoid the excitation of the terrestrial planets' eccentricities. 

This contingency requires that the delivery of LHB impactors to the Moon be achieved after the 
formation, large-scale migration, and/or global orbital restructuring of the giant planets. Within 
the context of the Nice model migration scenarios, only explaining the LHB requires that the bulk 
of giant migration occurred late. If giant planet migration was completed before the dynamical 
state of the inner solar system was determined, other aspects of the Nice model remain viable 
(e.g., the capture of Jupiter's Trojans and the orbital properties of the giant planets). Obviously, 
the delivery of the LHB impacto rs must then be accounted for via some other means (e.g., the 
Planet V hypothesis described in IChamberd 120071 ). 



Additional progress in unraveling the sequence of events and processes that produced the 
observed solar system might be made by examining the formation and early evolution of the gas 
and ice giant planets. We note that the nature and details of gas and ice giant formation in the 
solar system remain poorly understood. The basic assumptions and initial conditions of current 
migration models (e.g., that Jupiter and Saturn formed in a much more compact configuration) 
may change as the formation and early history of the giant planets are clarified. 

Finally, we note that migration near and through mean motion resonances is germane to the 
orbital evolution of planetary systems. Consequently, the resonant sweeping processes examined 
here have applications among extrasolar planetary systems, satellite systems, and other dynamical 
characteristics of the solar system. Analyses of these applications will be presented elsewhere. 
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